Fluctuations of local density of states and C speckle correlations are equal 
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We establish a conceptual relation between the fluctuations of the local density of states (LDOS) 
and intensity correlations in speckle patterns resulting from multiple scattering of waves in ran- 
dom media. We show that among known types of speckle correlations (C\, C'2, C3, and Co) only 
Co contributes to LDOS fluctuations in the infinite medium. We propose to exploit the equiva- 
lence of LDOS fluctuations and Co intensity correlation as a 'selection rule' for scattering processes 
contributing to Co. 



The local density of states of waves (LDOS) p(r, u>) 
is an important concept that regularly turns up in dis- 
cussions of waves in interaction with media. The num- 
ber p(r, uj)dVduj represents the local weight of all eigen- 
functions in the frequency interval du> around frequency 
u> inside a small volume dV around position r. In ho- 
mogeneous media it is independent of r and just equal 
to the 'density of states per unit volume' found in all 
textbooks. Near boundaries the LDOS exhibits Fricdel- 
type oscillations on the scale of the wavelength p|. In 
bandgap materials the LDOS was shown to govern the 
spontaneous emission of an atom at position r 0, 0] • In 
random media, where wave propagation is diffuse, the 
equipartition principle attributes an average local energy 
density of radiation that is directly proportional to the 
ensemble- averaged LDOS. This apparently simple princi- 
ple can have surprising consequences, for instance when 
waves with different velocities participate in the diffusion 
process, as is the case for seismic waves For disor- 
dered bandgap materials p| the equipartition principle is 
surprising in the sense that the multiple scattering pro- 
cess, with typical length scale the mean free path £ that 
is, in general, much larger than the wavelength or the 
lattice constant, distributes energy with sub-wavelength 
structure. From a fundamental point of view, the LDOS 
is also the crucial quantity in the recent studies on 'pas- 
sive imaging' Its basic principle is that for a homo- 
geneous distribution of sources — such as noise — the 
field correlation function (with time and space) is essen- 
tially proportional to the (Fourier transform of) LDOS, 
and thus sensitive to local structure, random or not. 

In random media the LDOS is a random quantity. Its 
statistical distribution has been studied previously within 
the frameworks of nonlinear sigma model 0, , random 
matrix theory 0, and optimal fluctuation method |ld| . 
The purpose of this paper it to establish a relation be- 
tween the fluctuations of LDOS — within the ensemble 
of random realizations — and intensity correlations in 
speckle patterns. Several contributions to intensity cor- 
relations have been identified. The 'standard', Gaussian 



correlation C\ is the best known ^lji but non-Gaussian 
correlations C2 and C3 have been predicted 0] an( i ob- 
served [13 , mostly in the transmitted flux. Recently the 
Co has been added 0E3. The C correlation is caused 
by scatterers close to cither receiver or source and is, 
surprisingly, of infinite spatial range. Contrary to the 
other correlations, Co is quite non-universal and highly 
dependent on details of the scatterers, such as their phase 
function. The total transmission coefficient is known to 
be dominated by C2, and the conductance by C3. Un- 
fortunately, the basic observable whose fluctuations are 
dominated by Co has never been identified. This is per- 
haps why its observation has never been reported so far. 

The fluctuations of LDOS can in principle be found 
from the average (Bethe-Salpeter) two-particle Green's 
function, but the diffusion approximation that is usually 
employed for this object 16] is not valid on length scales 
of the order of the wavelength, which appear to give an 
important contribution. Mirlin || noticed that in 3D 
the result is dominated by nearby scattering and is of 
order l/k£, where k is the wavenumber and £ 3> 1/fc is 
the mean free path. An exact calculation in the infinite 
medium with Gaussian white-noise disorder gives 
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Here G(r, r') is the Green's function of the wave equation 
describing the waves in the random medium. In the sec- 
ond equality we restricted to single scattering in the Born 
approximation [see Fig. [Ha)]. The value ir/k£ agrees ex- 
actly with the one found for the Gp speckle correlation 
[lower right diagram in Fig.^b)] [l4(. Going beyond sin- 
gle scattering, i.e. replacing the dotted lines in Fig. ^ by 
diffusion ladders, yields small corrections ~ l/(k£) 2 to 
both the variance of LDOS and Go- A deeper, generally 
valid relation between the fluctuations of LDOS and Go 
correlation is suggested by the above observations. This 
is the principal subject of the present work. 

Let us consider the simplest model possible: scalar 
waves in an infinite random medium with white-noise 
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FIG. 1: (a). Two-field intensity diagrams that give the lead- 
ing order to the variance of LDOS in a random medium with 
Gaussian uncorrelated disorder, (b). Typical four-field di- 
agrams that contribute to speckle correlations in a random 
medium. Solid and dashed lines denote retarded and ad- 
vanced Green's functions, respectively, shaded boxes are lad- 
der propagators (L) and Hikami boxes (H), dotted line with 
crosses denotes scattering of two wave fields on the same het- 
erogeneity. The variance of the local density of states (LDOS) 
can be obtained by integrating over x and x'. Only the lower 
right diagram yields a non- vanishing contribution to the vari- 
ance of LDOS in the infinite medium. 



disorder and leave more complicated situations for future 
work. Our assumptions are: 

1. At long distances, the diffusion approximation for 
the correlation (G(r, x)G*(r', x')) of two Green's 
functions is valid. 

2. The correlation (I(r, x)/(r', x')) of two intensities 
/(r, x) = |G(r,x)| 2 propagating from the source r 
to the receiver x is composed of terms of only four 
different classes, referred to as C\, G2, G3 and Go, 
distinguished by a different correlation range. 

The first assumption excludes 2D and ID random me- 
dia that are subject to localization effects. We will 
thus concentrate on 3D random media. The classifi- 
cation into Gi, G2, G3 and Go summarizes the out- 
come of numerous theoretical a ppr oaches and experi- 
ments [H IH 13 El 13 13 13H3 ■ The class Gi has 
short range correlation in both the source positions r, r' 
and the receiver positions x, x' (with a range at most 
equal to the mean free path). G2 has two parts. The 



first part has long range correlation (typically a power- 
law) in source positions and short range correlation in re- 
ceiver positions, the second part vice versa. G3 has long 
range correlation in both source and receiver positions. 
The class of terms described by Go exhibits an infinite 
ran ge c orrelation in either source or receiver positions 
[T3.ll5j . The classes G2, G3 and Go imply non-Gaussian 
statistics of the wave field. For weak disorder (kl ^> 1) 
this statistics is Gaussian and C\ dominates. 

The random dielectric constant is denoted by e(r), and 
we shall add a fictitious, homogeneous dissipation e a and 
later consider e a J. 0. The Green's operator for scalar 
waves is G(r, p, uj) = [(e(r) + ie a )ui 2 /c 2 — p 2 ] -1 - The re- 
solvent identity states that G — G* = —2ie a (uj 2 / c 2 ) GG* . 
In real space this translates to the identity 



ImG(r, r,uj,e a = 0) 



lim e a I dxl(r, x) (2) 



where the integral extends over the whole space, and 
the intensity J(r, x) was defined in assumption 2 above. 
We recall that the (radiation) LDOS p(r, uj) is equal to 
— (w/7rc 2 )ImG(r, r, uj) 0. Thus, Eq. expresses phys- 
ically that for a homogeneous distribution of sources, 
the local radiation density is directly proportional to the 
LDOS. For brevity we shall drop the frequency reference. 
The second moment of the LDOS can be expressed as 

(^W 2 ) = ^Ts^A / dx / dx'</(r,x)/(r,x')) (3) 

7T Z C° £a!0 J J 

Equation J2J establishes a conceptual relation between 
the variance of LDOS at r, Var[p(r)] = (p(r) 2 ) — (p(r)) 2 , 
given by its l.h.s., and intensity correlations in a speckle 
pattern created by a point source at r (the integrand of 
the r.h.s.). This facilitates a direct correspondence be- 
tween the various contributions to LDOS variance and 
speckle correlations. We demonstrate in Appendix 
that, among the four classes of speckle correlations, only 
Go contributes to the LDOS variance, because of it infi- 
nite range in the receiver positions x and x'. The con- 
tribution of Gi to Eq. © vanishes because Gi is short- 
ranged. The contributions of the long range correlations 
G2 and G3 vanish because they originate from crossings 
between diffuse propagators that respect current conser- 
vation [13. 

We conclude that the normalized fluctuations of LDOS 
and the Go speckle correlation are one and the same 
thing: 



Var[p(r)] 
<P(r)) 2 



= G 



(4) 



and that observational attempts to confirm the existence 
of Go should focus on the LDOS, either probed by spon- 
taneous emission or by using evanescent waves |l<Sj . 
It follows from our analysis that only correlations with 
infinite spatial range contribute to Var[p(r)], and Eq. (0J 
might serve as a definition for Gq. Alternatively, any 
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nonzero variance of LDOS implies the existence of spa- 
tial correlations of the intensity / with infinite range. 

Since Co correlation is non-universal and sensitive to 
the local, microscopic structure of the random medium, 
our Eq. Q implies that the fluctuations of LDOS are 
non-universal too, contrary to the universality of conduc- 
tance fluctuations. In the context of 'imaging with noise' 
essentially relying on the measurement of LDOS, the 
equivalence of Co correlation and LDOS fluctuations im- 
plies that only objects closer than a wavelength can affect 
the LDOS and can thus be imaged. 

The Co correlation determines the variance of LDOS 
at a given frequency to, and it continues to do so in the 
correlation of LDOS at two frequencies differing by some 
fl <C uj. We obtain (p(r, u)p(r, lu + fl)) c / (p(r, uj)) 2 ~ C , 
independent of Q. Similarly, if the disordered medium is 
not stationary, like, e.g., a suspension of small particles 
in Brownian motion, the LDOS will fluctuate in time. 
The time correlation of these fluctuations, (p(r, t)p(r,t + 
T)) c / (p( r, t)) 2 , is again determined by Co. According to 
Ref. [Tja], Co(t) decays as r~ 3 / 2 for large enough r. We 
conclude therefore that LDOS exhibits long-range cor- 
relations in time and infinite range correlations in fre- 
quency. 

To summarize, our main conclusion is that fluctua- 
tions of local density of states for waves in random media 
are conceptually equal to the recently predicted, though 
not yet observed Co intensity correlation, and not to the 
other known types of intensity correlation, which have all 
been observed. Crucial for this equivalence is the infinite 
spatial range of Co- In a finite medium the intensity cor- 
relations Ci.2,3 will emerge as extensive contributions to 
the LDOS variance, vanishing in some way as the medium 
scales upwards. Our analysis does not apply to localized 
media, where all correlation classes might contribute, in- 
tegrated over the finite volume £ 3 , with £ the localization 
length. With some minor modifications our main con- 
clusion should hold for infinite 3D disordered bandgap 
materials, where the LDOS is a much less trivial quan- 
tity. The equivalence between LDOS fluctuations and 
Co intensity correlations can serve as a selection rule for 
identifying scattering processes contributing to Co- 
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APPENDIX A 

In this Appendix we demonstrate that Ci, C2, and C3 
correlation functions do not contribute to the fluctuations 
of LDOS in Eq. (|3J), and that Co gives the only nonva- 
nishing contribution. We restrict ourselves to infinite, 
reciprocal media where G(x, r) = G(r, x) and assume 
H> I. 

We first consider Gaussian (Ci) statistics 



according to which (G(I)G*(2)G*(3)G(4)) 
(G(1)G*(2))(G*(3)G(4)) + (G(I)G*(3))(G*(2)G(4)). 
The first term just gives the average LDOS squared. 
In the diffusion approximation (assumption I), the 
correlation of two Green's functions takes the form [2l| : 

(G(r,x)G*(r',x')> = 

= (-Im G(r, r')> L(r, x) (-Im G(x, x')> (Al) 

In the infinite medium, the field propagator (ImG(r,r')) 
oscillates algebraically on the scale of the wavelength and 
decays exponentially beyond the extinction length i. The 
ladder propagator L(r,x), however, is very long-range 
and decays significantly only by absorption. Therefore, 
for the purpose of this paper we do not have to discrim- 
inate between r and r' or x and x' in L. On long length 
scales L obeys a diffusion equation with absorption time 

-L>V 2 L(r,x) + — L(r,x) = K6(r-x) (A2) 

where the factor K = lim ea ^o [ 7r (/ 5 ( x )) w £aii] 1 is im- 
posed by the ensemble-average of Eq. (0) 22J . 

The variance of LDOS caused by C\ becomes [see the 
upper left diagram in Fig. ^b)] 

Vari[p(r)] = ^ lune 2 a { p ( r )) 2 [ dxL 2 (r,x) 

x J dAx(-ImG(Ax)) 2 (A3) 

The integrand of the second integral is 
sin 2 (fcAa;) exp(— Ax/£)/(Ax) 2 , making the integral 
converge after typically the extinction length £, without 
the need for absorption. The integrand of the first in- 
tegral is typically L(x) 2 ~ (K 2 /D 2 x 2 ) exp(-2x / y/Dr^) . 
The critical contribution of Eq. (|A3|) comes from large 
x which justifies the diffusion approximation employed 
here. The first integral thus scales as ^fr^ . Since 
r a ~ l/e a we conclude that as e a | 0, the C\ contribu- 
tion to the variance of LDOS vanishes. All diagrams 
with short-range spatial correlations in both source and 
receiver positions have the same fate, in particular the 
diagram C[ In Fig.^b), which we discuss below. 

We now turn to C2, the first non-Gaussian contribution 
to the intensity correlation ^| • This is caused by a single 
crossing at an arbitrary point s in the medium [see the 
second and the third diagrams in the left column of Fig. 
IHb)], and is described mathematically by the 'Hikami 
box' vertex, with a scalar constant Ti. that needs not be 
specified here. Two very similar contributions exist that 
differ only in selection rules |19| . The first is short-range 
for x/x', and equals 

</(r,x)/(r,x'))c 2o = W7r 2 (p(r)) 2 (ImG(x,x')) 2 

x J dsL 2 (r,s)(Vi • V 2 )L(ri = s,x)L(r 2 = s,x) 

(A4) 
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According to Eq. we need the double integral 
e 2 JdxJ dx! of this object and let e a tend to zero. 
One integral converges again rapidly after one extinction 
length and is finite without absorption. We shall write 
J G?x'(ImG(x,x')) 2 = Vq and rearrange expression (|A4|) 
to 

Var 2a [p(r)] = ^HV (p(r)) 2 



x lime 2 J c?x J dsL 2 (r,s) |VL(s,x)| 2 
- V \ime 2 a J dx|VL(x)| 2 J dsL 2 {s) 



(A5) 

In the second step we conveniently made use of trans- 
lational symmetry of the infinite medium. We see that 
|Vi(x)| ~ 1/x 2 for large x, making the integral con- 
verge without the need of absorption. Its divergence for 
x < t is an artifact of the diffusion approximation which 
is not valid at small length scales, and which we shall 
ignore. Hence the volume integral over x is just finite, 
without absorption. The integral over the position s of 
the Hikami box scales as Jr^,. As e a | we conclude 
that the contribution of the first C 2 term to the LDOS, 
Var 2a [p(r)], vanishes. 

The second contribution from C 2 is long-range as a 
function of x — x' 20] . Its expression reads 

(7(r,x)/(r,x'))c 26 -^ 4 (p(r)> 2 (p(x)> 2 

x J ds L 2 (r, s)(Vi • V 2 )£(ri = s, x)L(r 2 = s, x') 

(A6) 

and a little rearranging shows that its contribution to the 
variance of LDOS is 



Var 2b [p(r)] 



lim e 2 



dxVL(x) 



dsL 2 ( 



(A7) 



We note that J v dxVL(x) = J A t v \ rfAL(x), where A(V) 
is the surface enclosing the volume V. The surface inte- 
gral vanishes for any closed surface since £(x) does not 
depend on the direction of x. Thus, Var 2 &[/9(r)] = 0- 

The contribution of C3-correlation, the origin of uni- 
versal conductance fluctuations, can be handled similarly. 
C3 contains two Hikami boxes [C3 in Fig.^b)], but that 
is a technical complication, and in just the same way as 
for C 2 it can be shown to vanish as e a J, 0. The diagram 
C'i in Fig. [Hb) looks very much like C3 but has actually 
short spatial range in both source and receiver positions. 
As a result it belongs to the class C\ , and its contribution 
to the variance of LDOS vanishes for the same reason as 
seen in Eq. (|A3|) . 

Finally, the Co correlation is given by [see the lower 
right diagram in Fig. ^b)] 

(7(r,x)I(r,x')) Co = (I(r, x)> (J(r, x')> x C (A8) 

with Co a dimensionless scalar depending on the nature 
of the scatterers. For weak white-noise, uncorrelated dis- 
order Co = n/k£ 0, 0. The essential property of Co 
that is important here is its infinite spatial range caused 
by the scattering of waves going to arbitrarily distant x 
and x' on a common scatterer in the vicinity of the source 
at r. Inserting Eq. (|A8|) into the expression for the LDOS 
variance and making use of Eq. IjAlfl we obtain 



Var [p(r)] = uj 2 : C q tt 2 (p(r))' 



x lim e 2 



dx (p(x)> L(r, x) 



C^ 2 (p{v)) 2 {p{x)Y lim e 2 a uj 2 K 2 T 2 
Co( P (r)) 2 



(A9) 



Hence Co provides the only surviving contribution to 
Var[p(r)]. 
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